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Abstract 

Perturbation expansions up to third order for the generalized spiked harmonic os¬ 
cillator Hamiltonians H = + x 2 + (A > 0, 2 y > a ,7 = 1 + \\/\ + 4 A), and 

small values of the coupling A > 0 , are developed. Upper and lower bounds for the 
eigenvalues are computed by means of the procedure of Burrows et al [J. Phys. A: 
Math. Gen. 20 , 889-897 (1987)] for assessing the accuracy of a truncated perturbation 
expansion. Closed-form sums for some related perturbation double infinite series then 
immediately follow as a result of this investigation. 
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1. Introduction 


It is well known that, although many perturbation expansions diverge, they may 
actually be asymptotic expansions whose first few terms can yield good approxima¬ 
tions. The family of spiked harmonic oscillator Hamiltonians 

d 2 X 

H = H 0 + XV = ——j + x 2 H-(0 < x < oo) (1.1) 

dx z x a 

affords interesting examples of this phenomenon. Harrell [1] have shown that the fa¬ 
miliar Rayleigh-Schrodinger perturbation series diverge according as n > where 
n is the order of the Rayleigh-Schrodinger term. For example, the first-order pertur¬ 
bation correction diverges for a > 3 , while the second-order correction term diverges 
if a>§, and so on. In a sequel of articles, Aguilera-Navarro et al [2], Estvez-Breton 
et al [3] , and Znojil [4] have shown for the case of a < 5/2, the so called ‘non-singular’ 
case, that the perturbation series of the ground-state energy up to the second-order 
corrections is given by 


E(X, a) = 3 + 
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for a < 5/2. (1.2) 


Based on resummation techniques, an analysis of Aguilera-Navarro et al [2] showed 
that 
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(1.3) 

where 2 Fi(a,b;c-,z) is the known Gauss hypergeometric function [5] with circle of 
convergence \z\ = 1. For the limiting case a —> 2 , the Erst term on the right-hand side 
of (1.3) was shown by Estvez-Breton et al [3] using l’Hopital’s rule to be 
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(1.4) 


Znojil, soon afterwards [4], showed elegantly that (1.4) follows immediately by ma¬ 
nipulating the Maclaurin expansion of the gamma function. Recently, Hall and Saad 
[6 - 10] investigated a larger class so called generalized spiked harmonic oscillator 
Hamiltonians 

d 2 . a \ 

(1.5) 


H-Ho + XV- —^+x 2 + ^ + — {A > 0). 


The Gol’dman and Krivchenkov Hamiltonian H 0 = --£? + x 2 + which admits the 
exact solutions 


4>n(x) = (- 1 )" 


2(7) 


n!r(y) 


^-i e -i x \F 1 (-n,'i,x 2 ), 


( 1 . 6 ) 
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with exact eigenenergies 


E n = 4 n + 2y, n = 0,1, 2,, 


7=1 + -Vl + 4A, 


(1.7) 


is regarded as the unperturbed part, and the operator V{x) = x~ a as the perturbed 
part. They obtained [8] the energy expansion up to the second-order as 


E( A, a) 
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r 2 ( 7 ) 
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( 1 . 8 ) 


valid for a < 7 + 1, where 7 = 1 + |\/1 + 4 A. A closed form sum for the infinite series 
in (1.2) appears as especial case. In particular, for 7 = 3/2 or A = 0, Eq.(1.8), for 
a < |, reduces to 


3 — a 
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and closed-form sums of the infinite series in (1.2) follow immediately. Furthermore, 
for a = 2, since 


4F 3 (1,1,2,2;2,2,7+1;1) = 2 F 1 (1,1;7+1;1) = 


r( 7 + i)r( 7 -i) 
r( 7 )r( 7 ) 


( 7 - 1 )’ 


by means of Chu-Vandermonde theorem [5] 


2+1 (a, b; c; 1 ) 


r(c)F(c — a — b) 
F(c — a)T(c — b ) ’ 


for c - a - b > 0, 


(1.9) 


the perturbation expansion (1.8) takes the very simple form 

E(X, a = 2) = 2 7 + —— - 4( ^'_ 1)3 + ■ ■ ■ (1-10) 

This is obtained, as expected, by means of Taylor’s expansion of the exact energy 
2 + ^/1 + 4 (A + A) about A = 0. In order to understand the result (1.4), however, we 
should note first 
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where we have used (1.9). Now since 
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we have 


y- (1)? 

^ (* + 1) (7)* * ! 
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(7 — 1)V’ (1) (7 — 1) — 1 


for 7 > 1, 


( 1 . 11 ) 


where ip^\z) is the first-derivative of the digamma function (or logarithmic derivative 
of the gamma function [ 11 ]) . Further, since V^XI) = tt > Hie result of (1.4) follows 
immediately by replacing 7 with 3/2 in (1.11). 

The interesting feature of the expression (1.8) is that, it can be applied to the 
ground-state eigenenergy at the bottom of each angular-momentum subspace labelled 
by l = 0,1,2,. .. in AT- dimensions: we just need to replace A with A —> A + (l + |(iV - 
1 ))(Z + i(AT - 3)). Furthermore, as we shall prove in the next section, for a = 4 and 
7 > 3 (or A > 3.75), the perturbation expansion (1.8) takes the very simple form 


E( X, a 


4) = 2 7 + 
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47 2 — 157 + 13 
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( 1 . 12 ) 


where ip is the digamma function. For a = 6 and 7 > 5 (or A > 15.75), (1.8) becomes 
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(1.13) 


In Sec. 2 we shall extended these perturbation expansions to third-order corrections. 
In Sec. 3, we shall discuss upper and lower bounds for the eigenvalues by means of 
the procedure of Burrows et al [ 12 ] for assessing the accuracy of a truncated pertur¬ 
bation expansion. These bounds will shed some light on the question regarding the 
acceleration of the variational method. Our conclusions and some remarks concerning 
the sums of some double infinite series will be given in Sec. 4. 

The functions 1 F 1 and 4 F 3 , mentioned above, are special cases of the generalized 
hypergeometric function [13] 
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where p and q are non-negative integers, and none of the fy, (j = 1 , 2 ,..., q) is 
equal to zero or to a negative integer. If the series does not terminate (that is to 
say, none of the a % . i = 1,2, ... ,p , is a negative integer), then the series, in the case 
p = q + 1 , converges or diverges accordingly as b| <1 or \z\ > 1 . For z = 1 , the 

q p 

series is convergent provided J2 ~ Here (a) n , the shifted factorial (or 

1=1 *=1 

Pochhammer symbol ), is defined by 

r(a + n) 

r(a) 


(a)0 = 1, (a) n = a(a + 1 )(a + 2)... (a + n - 1) 


n = 1,2,.... 


(1.15) 
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2. Third-order perturbation expansions 

In this section we will expand the perturbation expansions (1.8) to the third-order 
correction. Although, we will concentrate on the cases of a = 4 and a = 6 , since they 
are the most relevant in the literature [15 - 24], for other values of a the procedure 
is similar. In order to lay the foundation of the perturbation expansion (1.8), we first 
review the Rayleigh-Schrodinger perturbation theory for a non-degenerate case [25]. 
The fundamental problem in perturbation theory is the solution of the Schrodinger 
equation Hcj) = E( A)</> when H = H 0 + XV. The basic assumption is that b and E( A) 
may be expanded in power series in the perturbation parameter A: 

OO OO 

b = bo + E A *^’ E(X) =E 0 + ^2\ i e i . (2.1) 

i =1 i =1 

Here bo is a solution to the unperturbed problem H 0 ip 0 = E 0 ip 0 • We also choose the 
normalization (bo, b) = 1 , which implies that the higher-order corrections fa, fa... are 
orthogonal to bo • Perturbation theory tells us in this case that 

ei = (bo, Vfa), £2 = (bo, V fa), e 3 = {fa, V fa) — £\{fa, fa), ..., (2.2) 

or, equivalently [26], 


ei = (bo, Pbo), C2 


^ IHhI 2 

Ei — Efa £3 


EE 
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IPod 2 
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(2.3) 


From (2.2) it is clear that the first-order wave function fa determines the energy to 
the third-order. The matrix elements V l3 = (b*, Vifa) in (2.3) are computed by means 
of the basis solution {ip n } of the unperturbed Hamiltonian H 0 . For the generalized 
spiked harmonic oscillator Hamiltonian (1.5), the expectation values of the operator 
V(x) = x~ a with respect to the Gol’dman and Krivchenkov basis (1.6) are given 
explicitly by 


Vij = (- 1 ) 
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Of particular interest is 
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Recently, Hall et al [27 - 28] have shown that the first-order correction of the wave- 
function, in the case of a = 2 , is given by 
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Therefore from (2.2) and (2.3), by using (2.5), we have 
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as shown previously using summation technique. This idea can be used to obtain 
a simple form by expressing 4 F 3 in (1.8) in terms of elementary functions. These 
indeed are facilitated by the closed expression of the first-order correction of the wave 
functions developed earlier [27 - 28]. In the case a = 4, the first-order correction of 
the wave function reads 


M x ) 


1 x 1 2 e *2 

2v / 2 (7 — 2)(7 — 1) v^rby) 
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x 2 



for 7 > 2 . 


(2.7) 


where is the digamma function [ 11 ]. Using (2.2) and (2.3), we have 


4 F 3 (1,1,3,3;2,2,7+1;1) 


1 7(47 2 — 157 + 13) 
4(7- l)(7-2)(7 —3)’ 


7 > 3 


and therefore the perturbation expansion (1.12) follows immediately. These particular 
values of 4 F 3 ( 1 , 1 ,3,3; 2 , 2 , 7 + 1 ; 1 ) can be verified by means of the following lemma that 
extends the earlier identity 


3 F 2 (a,b,c + l;d, c;z) 


2-Pi(a, 6; d\ z ) H- - z 2 Fi(a + 1, b + 1; d + 1; z) 

cd 


given by Luke [29]. The proof follows immediately by use of the series representation 
for the hypergeometric functions 3 F 2 and 2 F 3 , as given by (1.14). 


Lemma 1: For M < 1 , 


4 E 3 (a, b,c+l,d+l; e,c,d; z) = 2 F 1 (a, 6; e; z) + — (1 + —)z 2 F 1 (a + l,b+l;e + l;z) 

ec d 

{a) 2 (b ) 2 2 


( 2 . 8 ) 


H- 7 7 \— 2F1 (cl + 2, b + 2; e + 2; z) 

dc(e )2 


Further, in the case of \z\ = 1 and e - a - b > 2 

T(e) 


4 J 7 3 (a, b, c + 1, d + 1; e, c, d; 1) = 


r(e — a)r(e — 6) 

(a) 2 (b ) 2 


r(e — a — b) 4-(1 H——)F(e — a — b — 1) 

c a 


+ l^g r (e- a -U 2 ) 
dc 


(2.9) 


In the case of a = 6 , the first-order correction of the wave function reads [27 - 28] 
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consequently, from e 2 = (if>o,x 6 fa ), we have for 7 > 5 


4+3 (1,1, 4, 4; 2, 2, 7 + 1; 1) 


7 / (7 — 2 )( 7 — 1) 2 ( 7 - 1) (40 — 57 7 + 24 7 2 — 3 7 3 ) \ 

18 \ (7 — 5) (7 — 4) ( 7 - 4 ) ( 7 - 3) ( 7 - 2) ( 7 - 1) ) 


( 2 . 11 ) 


Therfore Eq.(1.8) takes the simpler form (1.13), as the result of (2.11). In order to 
extend (1.12) and (1.13) to the third-order perturbation correction, we need only use 
the expression e 3 = (fa,Vfa) - ei(0i,0i), as mentioned in (2.2). Before we proceed 
with our calculations we shall first prove the following general result concerning the 
first-order correction of the wave function. 


Lemma 2: The first-order perturbation correction fa{x) °f the exact solution of 
Hamiltonian (1.5), with arbitrary a, satisfies the following normalization condition 


(0i) 0i) 


a 2 r 2 ( 7 -f) 
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as long as a < 7 + 2 . 

PROOF: We note that, by comparing the expression for e 3 in (2.2) and (2.3), we find 


( 01 , 0 l) 
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For the Hamiltonian (1.5), Vo » is given by (2.5) and E t is given by (1.7); therefore 
we have 


( 0 i, 0 i) 
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5 ^ 4 ( 1 , 1 , 1,1 + 1 ,| + 1 ; 2 , 2 , 2 , 7 +l;l) 


where we have used the Pochhammer identities (a) n+ i = a(a + l) n , (l) n = n\ and 
( 2 )„ = (n+ 1 )! (see (1.15)), and the series representation for the hypergeometric function 
5 F 4 , as given by (1.14). □ 


Direct computations, using e 3 = (</>i,a: _4 0i) -ei(0i,0i) where fa is given by (2.7) 
and ei = F p 7 ( ~ 2) leads, for a — 4 and 7 > 4 , to 


E{ A, a 


4) =2 7 + 


+ 


A 2 47 s — 157 +13 

(7-l)(7-2) 4(7 -1) 3 (7 -2) 3 (7 -3) 

I 67 5 - 1757 4 + 742 7 3 - 15257 2 + 1520 7 - 590 
8(7-4)(7-3) 2 (7-2)5( 7 -l) 5 


A 3 + ... 


( 2 . 12 ) 


For the case of a = 6 , the first-order correction of the wavefunction is given by (2.10). 
After some straightforward algebraic calculations, the ground-state perturbation ex¬ 
pansion, up to the third-order of A and valid for 7 > 7, now reads 


E{ A, ol — 6) — 27 + eiA + e 2 A 2 + e 3 A 3 + ..., 


(2.13) 
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where 

__A_ 

61 (7 — 1)(7 — 2) (7 — 3) ’ 

= r 2 (—3 + 7 ) / (7 — 2)( 7 — 1 ) 2 ( 7 - 1 ) (40 — 37(19 + (—8 + 7 ) 7 )) A 

8r 2 (7) U7-5)(7-4) ( 7 - 4 ) ( 7 - 3)( 7 - 2 ) ( 7 - 1) / ’ 

and £3 = 7 ;, for 

h = 192088 - 6559057 + 9458117 2 - 751923 7 3 + 36 O 8 II 7 4 - 1071517 5 + 19257 7 6 - 19177 7 + 8 I 7 8 , 
h = 8(7 - 7) (7 - 5 ) 2 ( 7 - 4 )(7 - 3 ) 5 ( 7 - 2) 5 ( 7 - l) 5 . 

A first reading of the articles by Sinanoglu [30] (the main results of which are 
not affected by his false claim), or even the work of Morse and Feshbach [31] on 
perturbation theory, one understands that the expressions (2.12) and (2.13) are upper 
bounds to the exact energy since all the odd-order energies would form upper bounds 
to the exact energy. This is not in fact true because £2 in the general perturbation 
expansion ( 2 . 1 ) will always have a negative sign, thus not guaraneeing the upper 
bounds [32-33]. However, it is possible to obtain a definite upper bound to the exact 
eigenvalue by means of the perturbation expansion. Thus 

E(X, a)=E 0 + ei A + 1 ’ ( 2 - 14 ) 

where (</>i,bi) is given by Lemma 2. The upper bound (2.14) can easily be demon¬ 
strated by applying the variational principle to the approximate wave function b = 
V’o + A^i, where bo and bi satisfies the zero-and first-order perturbation equations 

Hobo = E 0 tp 0 , (H 0 - F 0 )bi = (Fi - F)bo- (2.15) 

In Table (1), we compare the upper bounds obtained by means of (2.14) in the 
case of a = 4 and those of Aguilera-Navarro and Koo obtained by variational analysis 
using appropriate trial functions. In this next section, we shall obtain the symmetric 
lower and upper bound by means of the method of Burrows et al [ 12 ]. 

3. Lower and upper bounds 

It is natural to ask: how small A should be for the perturbation expansions (2.12) 
and (2.13) to be valid?. The question can be answered by studying upper and lower 
bounds to the eigenvalues. Based on the difference between the bounds we can infer 
a def ini te indication of the accuracy of truncated Rayleigh-Schrodinger perturbation 
series, such as (2.12) and (2.13). Wide bounds show that the truncated Raylcigh- 
Schrodinger perturbation series is suspect, while tight bounds demonstrate the high 
accuracy of the truncated expansion. For our purposes, the most suitable procedure 
developed for assessing the accuracy of a truncated perturbation expansion is due to 
Burrows et al [ 12 ]. A brief review of the method is presented here: for further details 
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the reader is referred to the original article. Most derivations of bounds for eigenvalues 
of self-adjoint operators start from a consideration of positive definite function 

(M0> e), p(fa e)) = ([H - e\fa [H - e]fa = Hfa - (0, Hfa 2 + (e - (0, Hfa) 2 > 0, (3.1) 

where H is the operator in question, e is a positive parameter, and 0 is a suitably 
chosen (normalized) function. If we expand the normalized function (j> in terms of 
the complete set of eigenfunctions {fa} of H with eigenvalues E n { A), </> = > 

a n = {4>,fa), fa = 1 = E„ °'i ) we can express the positive dehnite function in (3.1) 
as 

=^al(E n ( A) - e) 2 > 0 

n 

Let us assume that we have picked the value of e to lie closest to the value of the 
i th eigenvalue E t , i.e. 

(/r(0, e), n(<t>, e)) = ^ a l( E nW ~ > ( E iW - e) 2 > 0 (3.2) 

n 

Combining (2.6) and (2.7), it can easily be seen that 

/-(e) < Ei(X) < /+(e), (3.3a) 


where 

f± (e) = e ± vll ffb II 2 -(</>, ^) 2 + (e-(b,ffb)) 2 - (3-36) 

It is not hard to show that f±(e) is indeed a monotonic increasing function of e. This 
result will turn out to be useful in the following discussion. The bounds of Burrows 
et al follow [ 12 ] by setting 


H(4>,E P {\)) = [H 0 + XV - E p {\)](f) 


(3.4) 


where 

p 

0 = A^i(^o + Abi), E p {X)=Eo + ^2x i € i (3.5) 

i =1 

and, for all p < 3, tp 0 and fa satisfy the zero- and first-order equations of the 
Rayleigh-Schodinger perturbation theory (2.15). Further, the = 1,2,3 are given 
by means of (2.2). Here, IVi in Eq.(3.5) is a normalization constant for the truncated 
first-order expansion of the exact wavefunction: 

Ah = (i+A 2 (b!,bi)r 1/2 . 


If b and E P (X) were exact, n = 0 . Thus we expect n to be small if (j> and E p ( A) 
are good approximations to the exact solutions. Consequently, a good test of the 
approximations (3.3a-b) may be made by examining the value of the norm || p, ||= sfji? . 
Simple calculations, using (3.5) and ( 2 . 2 ), now give 

|| nifaE^X)) ||= Ah A 2 (0i, (V - ei) 2 0i) 1/2 , 


(3.6) 




Energy bounds for a class of singular potentials and some related series 


page 10 


|| m(< i>,E 2 (\] )) ||= {|| 7,(0,ERA)) || 2 ±A 4 e 2 {e 2 - 27V 2 (e 2 + Ae 3 )}} 1/2 , (3.7) 

and 

M0,£ 3 (A)) ||={|b(0,E 2 (A)) || 2 +A 4 {2A 3 e2 e 3 (l-fV 1 2 )-A 2 (l + A 2 )fV 1 2 ^ + A 4 e2}} 1 / 2 (3.8) 

where we have re-produced the formulas of Burrows et al [ 12 ] for computational 
convenience. In this case, (3.3a) implies 

/_(E P (A))<E(A)</+(E P (A)) (3.9a) 

where 

f±(E P (X)) = E P (A) ± || g(<j>,E p ( A)) || for p = 1,2,3. (3.96) 

The only new integral (beyond the usual integrals of Rayleigh-Schrodinger perturba¬ 
tion series) is seen to be ( 0 i|(V - ei) 2 (R) which restricts the value of 7, for example 
in case of a = 4, to be greater than 4 even if we have used the first-order approx¬ 
imation ei (for which 7 > 2 is sufficient). This is, of course, due to the bound’s 
dependence on e 3 which required 7 > 4. The result in this case, however, is very 
useful [23,34] when the radial Schrodinger equation is characterized by large angu¬ 
lar momenta l . For 7 = 4.5 ( i.e. A = 12 or l = 3 for A = 1(1 ± 1)) and A = 0.001, 
the first-order perturbation correction yields 9.000 114 285 with an error bound of 
±4.8346 x 10 -8 . The second-order perturbation corrections yields 9.000 114 279 with 
error bounds of ±4.7879 x 10 -8 ; while e 3 yields 9.000 114 279 with an upper bound of 
9.000 114 327 and a lower bound of 9.000 114 231. Now, for any fixed 0, the bound¬ 
ing functions f±(E p ( A)) are easily shown to be monotonic increasing functions of 
E p ( A), p = 1,2,3, as we indicated above. Consequently the optimal bound for the set 
{ER A) = E 0 + Xe 1 ,E 2 {X) = E 0 ± Aei ± A 2 e 2 , E 3 ( A) = E 0 + Xe ± ± A 2 e 2 ± A 3 e 3 } is indeed given, 
for A < -^ , by 

f-(Ei{X)) < E(X) </+(E 2 (A)). (3.10) 

The inequality A < ^ allows us to order the approximated eigenvalues as ERA) > 
E 3 (A) > ERA), for the sign of e 2 is always negative and the sign of e 3 is positive for 
moderate values of A. In Table II we have verified these results by obtaining upper 
and lower bounds for the eigenvalues by means of (3.6-8); underlined values are the 
optimal bounds. Similar bounds can be obtained for the case of a = 6 by using (2.13). 
Although, the upper bounds obtained by this method arc less accurate than the upper 
bounds obtained by means of (2.14), the advantage of this method is the symmetric 
lower and upper bounds available through (3.9). 

4. Conclusions and some remarks 

The main results of the present article are concrete upper- and lower-bound for¬ 
mulas (2.14), (3.9), and (3.10). There are many variational methods avaliable to solve 
the eigenvalue problem for the Hamiltonian (1.5), however they provide only upper 
bounds and usually no information is avaliable concerning the accuracy of the method 
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other than comparison with numerical solutions of the Schrodinger equation in ques¬ 
tion. Furthermore, for very small values of the parameter A, variational methods are 
usually slow and a large number of the matrix elements are needed to obtain suffi¬ 
cient accuracy. We have presented upper and lower bounds for such situations which, 
as table (I) and (II) indicate, provide excellent results for very small values of A. 
Although, the techniques used to produce the present results are standard, the abil¬ 
ity of these techniques to generate explicit bounds is a consequence of our pervious 
achievements, yielding concrete forms for the first-order perturbation corrections of 
the wave functions. 

Aside from the upper and lower bounds obtained, there are also some interesting 
results concerning a closed-form sums for double infinite series that follow directly 
from the present work. It is clear from (2.2) and (2.3) that 


EE 


tbnhnmbmO 
(Eq ~ E n )(Eo — E m ) 


(biWbi) 


(4.1) 


where V nrn ,n = 1 , 2 ,...,m = 1 , 2 ,... are given by (2.4). We will now look at the cases 

a = 2 ,4, 6 ,_ Similar results can be obtained for the cases of a = 1 ,3,5,... by means 

of the first-order corrections for the wave functions given previously [27-28]; however, 
the calculations will be more involved for such cases. For a = 2 , we know that the 
matrix elements (i.e from (2.4)) read 


V n . 


t i \n+m C(7~1) / 

i L > r( 7 ) V "!( 7 )m 


if m > n, 


(-i rwa if n > to. 


(4.2) 


On other hand, the first-order correction of the wave function in this case reads 
[27 - 28] 


cj>i(x) = 


1 x 1 2 e 2 


iogOO - 2^(7) 


x/2( 7 -l)VT(7) 

Consequently, the following results follow immediately, 

Lemma 3. For 7 > 1 and V nrn as given by (4-2), we have 


for 7 > 1. 


(4.3) 


OO OO 


EE 


to n*nm hnO 
16nm 


+ 


V’ (1) ( 7 ) 


8(7 — l ) 5 16(7 — l ) 3 ’ 


(4.4) 


where ( 7 ) is the first derivative of the digamma functions. 

The proof of this Lemma is obtained by calculating the inner product of the right- 
hand side of (4.1) by means of (4.3) for 0 < x < 00 , where V(x) = x~ 2 and E n = 4n + 2y 
( n = 0,1,2 ,...). For the case a = 4 and 7 > 2 , the matrix elements (2.4) read 


f (—' 1 )" +m wyq I fjV // 7 $vfe[ 7 ( m -«+!) + 2n] if m > n, 

V nm =l / V E_ (4.5) 

[ (_l)n+mr| 2 z 2 |^/^^| 7 ( n _ m + 1 ) + 2 to] if n > to. 
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On the other hand, the first-order corrections of the wave function for this case are 
given by (2.7). Therefore (4.1) leads to the following results 

Lemma 4. For 7 > 4 and V nrn as given by (4-5), we have 


EE 


bo n bnm EnO 

16 nm 


-820 + 19547 - 17537 2 + 694 7 3 - 9O7 4 - 12 7 5 + 37® 
16(7 — 4) (7 — 3) 2 (7 — 2) 5 (7 — l) 5 


where 7 ) is the first derivative of the digamma functions. 


V ,(1) (7) 

16(7 — 2) 3 (7 — l) 3 ’ 
(4.6) 


As final case that we illustrate, namely a = 6 and 7 > 3 , we point to the fact that 
Eq.(2.4) lets us deduce 


t pn+m r(7-3) / m!(7)n N , 

1 2r( 7 +2) V n!(V)m X 

[(2 + m)(l + m ) 7(7 + 1) — 2n(l + m)( 7 — 3) (7 + 1) — n(l — n)( 7 — 2) (7 — 3)] if m > n , 


(- 1 ) 


n+m T(7-3) 




2r(7+2) Y m!(7)n 

l [(2 + n)(l + n) 7(7 + 1) — 2m(l + n )(7 — 3) (7 + 1) — m(l — m)( 7 — 2) (7 — 3)] if n > m. 

(4.7) 

where the Erst order correction for the wave function is now given by (2.10). Therefore, 
by means of (4.1), we conclude 


Lemma 5. For 7 > 7 and V nrn as given by (4-7), we have 

bonbnmLnO F ^ 


EE 

m =1 n= 1 


16n?n 


I 2 16(7 — 3) 3 (7 — 2) 3 (7 — l) 3 ’ 


(4.8) 


where 


h = 522652-17174407+237193l7 2 -17850467 3 +79206l7 4 -2069647 5 +287257 6 -11587 7 -1697 8 +167 9 , 
and 

h = 32(7 - 7)(7 - 5) 2 ( 7 - 4)(7 - 3) 5 ( 7 - 2) 5 ( 7 - l) 5 , 
where ^ ( 1 ) ( 7 ) the first derivative of the digamma function. 
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Table (I) A comparison between the upper bounds for the Hamiltonian (1.5), for 
a wide range of values of A = 1(1 + 1) and A, by formula (2.14) and the bounds E% 
obtained by Aguilera-Navaro et al [21]. Exact results E found by direct numerical 
solution of Schodinger’s equation are also presented. 


A 

l 

K 

E u 

E 

0.001 

3 

9.000 114 279 82 

9.000 114 279 12 

9.000 114 279 12 


4 

11.000 063 490 8 

11.000 063 490 7 

11.000 063 490 74 


5 

13.000 040 403 7 

13.000 040 403 6 

13.000 040 403 64 

0.01 

3 

9.001 142 268 25 

9.001 142 199 48 

9.001 142 199 40 


4 

11.000 634 795 5 

11.000 634 788 8 

11.000 634 788 89 


5 

13.000 404 001 8 

13.000 404 000 6 

13.000 404 000 60 

0.1 

3 

9.011 370 328 09 

9.011 364 261 69 

9.011 364 026 18 


4 

11.006 336 739 4 

11.006 336 100 1 

11.006 336 099 23 


5 

13.004 036 546 4 

13.004 036 432 5 

13.004 036 432 52 

1 

3 

9.109 013 250 38 

9.109 311 262 10 

9.108 658 607 52 


4 

11.062 293 143 4 

11.062 249 282 0 

11.062 241 719 38 


5 

13.040 025 483 8 

13.040 015 551 5 

13.040 015 183 06 


50 

103.000 400 037 

103.000 400 036 

103.000 400 036 76 


Table (II) 

_AL+z 2 +±4 


(3.10). 


Eigenvalue Bounds for different values of A for the Hamiltonian H = 
2* • The underlined values are the optimal bounds according to inequality 



£l 

£2 

£3 

A 

E l 

E u 

E l 

E u 

E l 

E u 

0.001 

9.000114234 

9.000114334 

9.000114231 

9.000114327 

9.000114231 

9.000114327 

0.01 

9.001138022 

9.001147691 

9.001137408 

9.001146987 

9.001137409 

9.001146989 

0.1 

9.010945111 

9.011912031 

9.010883476 

9.011841809 

9.010885097 

9.011843425 

1 

9.065963521 

9.162607906 

9.059599522 

9.155786288 

9.061245282 

9.157377241 
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